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Abstract 

Lg  spectra  from  20  imdergroimd  nuclear  explosions  and  52  shallow 
earthquakes  in  central  Eurasia  have  been  collected  at  21  broad  band  IRIS, 
CDSN  and  KNET  stations.  Using  the  non-hnear  method  of  Xie  (1993),  we 
simultaneously  invert  for  Lg  source  spectral  parameters  and  path-variable 
Lg  Qo  and  77  values  (Lg  Q  at  1  Hz  and  its  power-law  frequency  dependence, 
respectively).  The  inversions  yield  Lg  seismic  moments  (Mq),  comer  frequen¬ 
cies  (fc)  for  the  events,  as  well  as  Lg  Qo  and  77  values  for  numerous  paths  in 
central  Eurasia. 

Grossly  speaking,  Lg  Q  and  Lg  coda  Q  are  similar  in  central  Eurasia, 
although  minor  discrepancies  exist  between  (a)  Lg  77  and  Lg  coda  77  at  large 
(>  about  2700  km)  distances,  and  (b)  Lg  Qo  and  Lg  coda  Qo  values  in  a  sub- 
region  northeast  of  the  Lop  Nor  test  site.  For  both  explosions  and  earth¬ 
quakes  Lg  Mo  values  correlate  linearly  with  the  ISC  Mb  values,  both  having 
slopes  slightly  greater  than  1.0.  For  the  same  Lg  Mo  values,  Mb  values  from 
explosions  tend  to  be  larger  than  those  from  earthquakes.  Lg  Mo  tend  to 
scale  with  with  a  closer  to  4  than  to  3.  Regression  analysis  over  Mo  and 
fc  values  suggest  that  for  the  same  Mo  values,  explosions  tend  to  have 
higher  values.  This  may  form  a  basis  of  an  explosion  discriminant,  but 
the  fact  the  there  is  a  shght  overlap  among  Mo  and  fc  values  from  explosions 
and  those  from  earthquakes  suggests  that  the  use  of  this  discriminant 
should  be  used  with  caution. 

Applications  of  the  methodology  used  in  this  study  to  other  regional 
phases,  such  as  Pn,  may  contribute  to  our  understanding  of  regional  wave 
exdtation/propagation  by  various  types  of  seismic  source,  and  to  the  evalua¬ 
tion  of  various  discriminants  using  regional  waves. 
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Research  Accomplished 


Over  the  past  two  years  we  have  collected  Lg  spectra  from  20  under¬ 
ground  nuclear  explosions  in  the  Balapan  and  Lop  Nor  test  sites,  recorded 
by  17  broad-band  IRIS,  CDSN  and  KNET  stations  (Figures  1,  4,  also  cf.  Xie 
et  ah,  1995).  We  have  also  collected  Lg  spectra  from  52  shallow  earthquakes 
(5-33  km)  that  occurred  in  the  areas  of  the  central  Asian  Republics  of  the 
F.S.U.  and  southern  Xinjiang,  China,  recorded  by  11  broad  band  stations 
(Figures  1,  5).  These  spectra  are  used  to  invert  for  Lg  source  spectral 
parameters  (Mq,  fc)  and  path  variable  Lg  Qo  and  77  values  (Lg  Q  at  1  Hz  and 
its  power-law  frequency  dependence,  respectively),  with  the  non-linear 
inverse  method  of  Xie  (1993).  Figures  2  and  3  show  examples  of  the  fit  of 
the  optimal  source/path  parameters  to  the  observed  Lg  spectra,  where  the 
inverted  Qo,  77  values  for  multiple  stations  that  recorded  an  explosion  (Fig¬ 
ure  2)  and  earthquake  (Figure  3)  are  used  to  remove  path  effects,  resulting 
in  reduced  Lg  spectra  at  the  source.  It  appears  that  the  fit  for  the  explosion 
event  (Figure  2)  is  better  than  for  the  earthquake  event  (Figure  3),  particu¬ 
larly  for  the  station  averages  Gower  right  panels).  The  main  findings  of  this 
study  are: 

(1)  Importance  of  the  number  of  stations  recording  Lg:  During  the  spectral 
inversions  using  Lg,  we  foimd  that  when  the  number  of  stations 
recording  the  same  event  is  less  than  3,  the  available  Lg  spectra  is  typ¬ 
ically  not  sufficient  for  a  simultaneous  inversion  of  source  Mo,  fc  and 
path  Qo,  rj  values.  Accordingly,  a  priori  information  on  the  Qo,  77  val¬ 
ues  obtained  in  previous  simidtaneous  inversions  must  be  used.  It 
also  appears  that  the  Mq  and  fg  values  obtained  with  Lg  spectra  from 
only  one  or  two  stations  are  less  reliable  than  those  obtained  with  Lg 
spectra  from  more  stations.  Much  of  the  scatter  in  Figure  9  is  due  to 
the  Mo,  fc  values  obtained  iising  only  1  or  2  stations. 

(2)  Effects  of  radiation  pattern  by  earthquake  sources:  For  paths  connect¬ 
ing  the  Lop  Nor  test  site  and  stations  AAK,  GAR  and  ARU,  there  are 
two  sets  of  Lg  Qo  values  obtained,  one  using  earthquake  data  and  the 
other  using  explosion  data  (Figures  5  an  6).  These  Lg  Qo  values  are 
very  similar,  suggesting  that  the  radiation  patterns  for  the  earthquake 
sources  are  insignificant. 

(3)  Comparison  between  Lg  Qo  and  and  Lg  coda  Qo:  For  most  of  the  paths 
used  in  this  study,  the  Lg  Qo  obtained  in  this  study  are  highly  consis¬ 
tent  with  the  Lg  coda  Qo  map  of  Xie  &  Mitchell  (1991)  and  Pan  et  al. 
(1992)  (see  Figures  4,  5,  6).  For  three  paths  that  run  from  the  earth¬ 
quake  in  Xinjiang  northeastward  to  stations  MDJ,  HIA  and  TLY  (Fig¬ 
ure  6(a)),  the  Lg  Qo  values  are  somewhat  higher  than  the  Lg  coda  Qo 
values  (Figure  4).  The  most  hkely  cause  of  this  discrepancy  is  that  the 
direct  Lg  phase  and  Lg  coda  are  ciffected  by  3D  structural  complexities 
in  different  manners  in  the  area,  causing  the  two  Qo  values  differ. 
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(4)  Comparison  between  Lg  77  and  and  Lg  coda  tj:  When  the  epicentral  dis¬ 
tance  (A)  is  less  than  about  2700  km,  the  frequency  dependence  of  Lg 
Q,  T],  obtained  in  this  study  agrees  (within  an  imcertainty  level  of 
about  0.1  to  0.2)  with  the  Lg  coda  tj.  At  larger  distances  (A  >  2700  km), 
the  TJ  values  obtained  in  this  study  tend  to  be  low  (often  down  to  ~  0.0). 
This  discrepancy  is  most  likely  due  to  imprecisely  estimated  Lg  7  in 
this  study  due  to  narrower  frequency  bands,  or  effects  of  the  earth’s 
curvature  at  large  A. 

(5)  Scaling  of  Lg  Mq  with  ISC  Mj,:  For  both  explosions  emd  earthquakes, 
the  Lg  Mo  values  correlate  linearly  with  ISC  Mj,  (Figure  7).  Linear 
regression  over  the  points  in  Figure  7  3deld 

log  Mo  =  1. 19(±0. 11)  Mb  -t-  8. 85(±0. 64)  (1) 

for  explosions,  and 

log  Mo  =  1. 04(±0. 09)  Mb  -H  10. 66(±0. 51)  (2) 

for  earthquakes.  These  are  straight  lines  that  parallel  each  other  (Fig- 
me  7),  but  are  offset  such  that  for  the  same  Mo  values,  Mb  values  tend 
to  be  systematically  higher  for  explosions. 

(6)  Scaling  of  Lg  Mo  with  fc  for  explosions:  Figure  8  shows  Lg  Mo  versus  fc 
obtained  in  this  study  obtained  for  the  20  imdergroimd  nuclear  explo¬ 
sions  studied,  with  the  explosion  source  model  [i.e.,,  the  model  with  an 
overshoot  effect;  see  equation  (2)  of  Sereno  et  al.  (1988)  or  equation 
(10)  of  Xie,  (1993)].  Lg  Mo  values  correlate  linearly  with  ISC  Mb  and  a 
linear  regression  yields 

logMo  =  15.12(±0.22)-3.98(±0.43)logfc  .  (3) 

This  suggests  that  Lg  Mo  scales  with  instead  of  with  (i.e.,  con¬ 
stant  stress  drop  scaling). 

(7)  Scaling  ofLg  Mq  with  fc  for  earthquakes:  Figure  9  shows  Lg  Mo  versus 
fc  obtained  for  the  53  earthquakes  obtained  in  this  study,  using  an 
earthquake  source  model  (i.e.,  the  model  without  overshoot).  A  lin¬ 
ear  regression  over  the  points  in  Figure  9  yields 

log  Mo  =  14. 85(±0. 29)  -  3. 56(±0. 29)logfc  .  (4) 

(8)  Difference  between  the  Mo,  fc  scaling  for  earthquakes  and  explosions: 
Equation  (3)  is  obtained  for  Lg  Mo  and  fc  values  of  explosion  sources, 
obtained  in  inversions  using  the  explosion  source  model.  On  the  other 
hand,  equation  (4)  is  obtained  for  Lg  Mo  and  fc  values  of  earthquakes 
using  the  earthquake  source  model.  For  the  purpose  of  discriminating 
explosions  from  earthquakes,  it  is  desirable  to  obtain  Mo  and  fc  for 
explosions  using  the  earthquake  source  model,  thus  simulating  a  situ¬ 
ation  where  we  do  not  know  that  the  explosions  imder  study  are  explo¬ 
sions,  The  resulting  Lg  Mo  and  fc  values  for  the  explosions,  obtained 
using  the  earthquake  source  model,  are  plotted  in  Figure  9,  and  a  lin¬ 
ear  regression  over  these  values  yields 
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log  Mo  =  15.69(±0.23)-3.83(±0.45)logfe  .  (5) 

The  straight  hnes  represented  by  equations  (4)  and  (5)  are  subparallel, 
both  being  closer  to  Mq  ~  scaling  than  to  Mq  ~  scaling.  However, 
the  two  lines  are  offset  and  for  a  given  Mq,  the  explosions  tend  to  have 
higher  values.  This  suggests  greater  high-frequency  content  of  Lg 
from  explosions,  as  compared  to  that  from  earthquakes  of  similar 
moments.  The  Mq  and  fc  values  may  therefore  be  used  to  discriminate 
explosions  from  earthquakes.  However,  Figure  9  shows  that  there  is 
some  slight  overlap  of  the  Mq  and  fc  values  for  the  two  groups  of  events 
at  smaller  moments,  indicating  that  this  discriminant  should  be  used 
with  caution. 


Conclusions  and  Recommendations 

Grossly  speaking,  Lg  Q  and  Lg  coda  Q  are  similar  in  central  Eurasia, 
although  minor  discrepancies  exist  between  (a)  Lg  77  and  Lg  coda  77  at  large 
(>  about  2700  km)  distances,  and  (b)  Lg  Qo  and  Lg  coda  Qo  values  in  a  sub- 
region  northeast  of  the  Lop  Nor  test  site.  For  both  earthquake  and  explo¬ 
sion  sources,  Lg  Mq  values  correlate  linearly  with  the  ISC  Mb  values,  both 
having  slopes  that  are  slightly  greater  than  1.0.  Lg  Mo  tends  to  scale  with 
f^“,  with  a  being  closer  to  4  than  to  3.  Regression  analysis  over  Mq  and  fc 
values  suggests  that  for  the  same  Mo  values,  explosions  tend  to  have  higher 
fc  values.  This  may  form  a  basis  of  a  explosion  discriminant,  but  the  fact 
the  there  is  a  slight  overlap  among  Mo  and  4  values  from  explosions  and 
those  from  earthquakes  suggests  that  the  use  of  this  discriminant  should  be 
used  with  caution. 

Future  research  is  recommended  in  the  following  areas: 

(1)  Establish  more  precise,  perhaps  distance  and  frequency  dependent 
geometrical  spreading  for  the  Lg  phase  and  other  regional  phases 
based  on  synthetics  \ising  reahstic  velocity  structures. 

(2)  Conduct  experiments  to  see  if  the  co^  source  model  needs  to  be  modified 
for  Lg  excitation  by  earthquake  sources. 

(3)  Apply  the  same  methodology  in  this  study  to  the  spectral  characteris¬ 
tics  of  excitation  and  propagation  of  other  regional  phases,  particularly 
the  Pn  phase,  and  systematically  evaluate  the  Pn/Lg  spectral  ratio  dis¬ 
criminant.  Also,  test  to  see  if  the  Pn/Lg  discriminant  is  more  reliable 
than  the  Lg  discriminant. 
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Fig.  1.  Locations  of  the  20  underground  nuclear  explosions  (crosses),  52  earthquakes  (circles)  and 
21  seismic  stations  (triangles)  used  in  this  study.  The  numbers  of  stations  providing  Lg  records  are 
17  for  the  explosions,  and  11  for  the  earthquakes. 
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Fig.  2.  Synthetic  Lg  source  spectra  for  seven  KNET  and  IRIS  stations  recording  the  October  5, 
1993,  Lop  Nor  explosion,  versus  the  observed  Lg  spectra  that  are  reduced  to  source  by  removing 
path  effects.  The  lower  right  panel  is  the  average  for  all  of  the  seven  stations.  The  synthetic  spec¬ 
tra  are  calculated  using  optimal  source  spectral  parameters  (M©  =  8.3x10^®  Nm,  f;=0.68  Hz) 
obtained  in  the  inversion.  Path  Qo  values  obtained  in  the  inversion  are  written  on  the  top  of  the 
panels,  together  with  the  epicentral  distances. 


ARU,  6564.-147,  3006. 5km  AAK,  424  +  .93,  1301.3km  BJI,  393  +  -37,  2377.6km  STATION  AVERAGE 
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Fig.  3.  S3rnthetic  Lg  source  spectra  for  six  IRIS  and  CDSN  stations  recording  the  October  2,  1993, 
southern  Xinjiang  earthquake  (Mb  =  5.6),  versus  the  observed.  The  synthetic  spectra  are  calculated 
using  optimal  source  spectral  parameters  (Mo=2. 6  x  10^®  Nm,  f.=0.41  Hz)  obtained  in  the  inversion. 
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Fig.  4.  The  Lg  paths  for  which  Q  values  are  measured  in  this  study,  plotted  with  the  Lg  coda  Q© 
values  from  the  tomographic  inversion  by  Xie  and  Mitchell  (1991).  Solid  dots  and  triangles  are 
sources  and  stations,  respectively. 
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Fig.  5  Lg  Qo  values  obtained  for  the  great  circle  paths  from  the  Lop  Nor  and  Balapan  test  sites  to 
the  17  IRIS,  CDSN  and  KNET  stations.  Water-covered  areas  are  shaded. 


Fig.  6(a)  Lg  Qo  values 
and  CDSN  stations. 
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Fig.  6(b)  Lg  Qo  values  obtained  fo 
Republics  of  the  F.S.U.  to  IRIS  and ' 


mb  versus  Mo 


Mo  (10'*  Nm) 

Fig.  7.  Mb  values  versus  logarithm  of  Mo  values  (in  10^®  Nm)  obtained  for 
explosions  and  earthquakes.  Straight  lines  represent  the  linear  regression 
fitting. 


fc  versus  Mo  (Explosion  Model) 
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Fig.  8.  Logarithm  of  Mq  (in  10^®  Nm)  versus  logarithm  of  f<;  values  for  the 
explosions  studied,  obtained  by  inverting  the  Lg  spectra  using  the  explosion 
source  model.  Straight  line  represents  the  linear  regression  fitting. 
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fc  versus  Mo  (Earthquake  Model) 


Mo  (10**  Nm) 

Fig.  9.  Logarithm  of  Mq  (in  10^®  Nm)  verst:is  logarithm  of  f^  values  for  the 
earthquakes  and  explosions  studied,  both  obtained  by  inverting  the  Lg  spec¬ 
tra  using  the  earthquake  source  model.  Straight  lines  represent  the  linear 
regression  fitting. 
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